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Abstract. Visceral leishmaniasis (VL) seroprevalence in Kenya is unknown because of the lack of a practical and 
accurate diagnostic test or surveillance system. A nove! serological assay was used to estimate the seroprevalence of 
Leishmania-specific antibodies, and Global Information System and spatial clustering techniques were applied to study 
the presence of spatial clusters in Parkarin and Loboi villages in B aringo Districi in 2001. VL seroprevalences were 
52.5% in Parkarin and 16.9% in Loboi. Significant associations among seropositivity and house construction, age, and 
proximity to domestic animai enclosures were found. A significant spatial cluster of VL was found in Loboi. The spatial 
distribution of cases in the two villages was different with respect to risk factors, such as presence of domestic animals. 
This study suggests that disease control efforts could be focused on elimination of sand fly habitat, placement of domestic 
animai enclosures, and targeted use of insecticides. 



INTRODUCTION 

The leishmaniases are a group of mainly zoonotic infections 
caused by protozoan parasites of the genus Leishmania. 
These infections produce a variety of clinical diseases depend- 
ing on the virulence or tropism of the parasite and differential 
host immune responses.^ The World Health Organization es- 
timates that with 350 million people at risk worldwide, 12 
million people are infected with Leishmania parasites and 
that as many as 2 million new cases occur each year. Recent 
reports detail a visceral leishmaniasis (VL) pandemie in the 
Horn of Africa and parts of India, Nepal, Bangladesh, and 
Central and South America. A major epidemie reported in 
the Sudan from 1989 to 1993 was responsible for the deaths of 
approximately 10% of that country's population.^ 

Whereas the definitive diagnosis of VL leishmaniasis re- 
quires demonstration of parasites from smears or in vitro cul- 
tivation, these methods can be time consuming and involve at 
least one invasive, and at times risky, procedure such as 
lymph-node biopsy, splenic or bone marrow aspiration. A 
variety of serological methods, including indirect immuno- 
fluorescence (IFA),^ enzyme-linked immunosorbent assay 
(ELISA),"*'^ and polymerase chain reaction (PCR)'' are asso- 
ciated with a number of problems including cross-reactivity 
with other pathogens, high cost and/or need for sophisticated 
equipment. The direct agglutination test (DAT)^ is a sensi- 
tive, specific, and simple test, but the main disadvantage is 
that it can have high background seropositivity in endemie 
areas. 

New diagnostic tests for VL have been a major focus of 
numerous research groups.^"' Martin and co-workers re- 
ported the use of a soluble antigen prepared from Leishmania 
dono vani as the foundation for an ELISA that could detect 
specific IgG antibodies in kala azar patients.^° Soluble anti- 
gens produced by Leishmania (exo-antigens) are released 
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info the vertebrate host and sand fly vector. The assay was 
improved and developed info an ELISA that can detect 
Leishmania-specific IgM and IgG antibodies in patients with 
visceral and cutaneous leishmaniasis. '^^ In preliminary studies, 
the sera from 129 visceral leishmaniasis patients (Brazil, Italy, 
North Africa, Nepal) with sera from matched controls were 
tested. The test reported a sensitivity of 98.2%. Testing this 
ELISA with small subsets of naìve North American normals, 
Kenyan endemie normals, malaria, African trypanosomiasis, 
echinococcosis, filariasis, and schistosomiasis case patients 
yielded a specificity of better than 97%.^°"^^ 

Visceral leishmaniasis, caused by Leishmania donovani, is 
endemie in B aringo Districi, Kenya. The first case of VL was 
recorded in 1948 in the Districi Hospital at Kabarnet.^'' Some 
scientists believe that nomadic Turkanas may have intro- 
duced the disease info the area from the north. Others specu- 
late that Kenyan soldiers returning from North Africa after 
Word War II were responsible for introduction of the para- 
site. Phlebotomus martini is the vector of the parasite, and 
man is the only known reservoir.^"* Currently, the disease is 
widespread throughout Baringo Districi. 

A number of studies have used diagnostics, epidemiologie 
surveys, and Global Information System (GIS) technology 
to identify potential risk factors with several tropical dis- 
eases. ^^"^^ The current study applied a multidisciphnary ap- 
proach with integration of new technologies. The objectives 
in this study were to determine in two villages in a VL- 
endemic area of the Baringo Districi of Kenya 1) the current 
seroprevalence of Leishmania-specific antibodies, 2) the pres- 
ence of spatial clustering (hot spots) of seropositive cases, and 
3) the Identification of associations with certain risk factors in 
those same villages. 



MATERIALS AND METHODS 

Study area. Baringo Districi is the only region in Kenya 
where both VL and cutaneous leishmaniasis (CL) have been 
found together.^® It covers an area of approximately 10,000 
km^ and is located north of the Equator in Kenya's Rift Val- 
ley Province. Marigat Township is centrai to this region and a 
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site where most VL studies bave been conducted in the past. 
Generally, the district includes a range of landscapes varying 
between fertile highlands as high as 2,700 m elevation and 
semiarid lowland at about 900 m elevation. The Tugen hills 
and the adj acent highlands in the southwest receive 1,200- 
1,500 mm average yearly rainfall and bave daily average tem- 
peratures ranging from 10°C to 32°C. In contrast, lowland 
areas to the west, east, and north of the Tugen Hills receive 
annual rainfall from 300 to 700 mm and the temperatures vary 
between 16°C and 42°C. The rainy season is from March to 
September, with maximum rainfall in May and August. Three 
main ethnic groups, ali classified as agropastoralists, inhabit 
Baringo District: the Tugen, Pokot, and Njemps. The district 
is sparsely populated due to harsh physical and climatic con- 
ditions. There have been no reports of cutaneous or mucocu- 
taneous leishmaniasis in the two study areas. 

Study sites. This epidemiologie study was conducted in the 
villages of Parkarin and Loboi (0°35'N latitude, 36°06'E lon- 
gitude) in the Rift Valley Province, Baringo District of centrai 
Kenya, and approximately 207 kilometres northwest of Nai- 
robi. The population of Parkarin (A'^ = 286 inhabitants) was 
spread over an area of 3.3 km^. The population of Loboi {N 
= 643 inhabitants) was spread over an area of 2.3 km^. Both 
villages, being only 7 km apart, have the same subtropical 
climate. Outwardly, the two villages appear to be very similar 
and both have had a long history of kala azar. The villagers 
live in family units (homesteads) each consisting of several 
huts. Within the homestead, householders may Uve together 
in one house or in 2-3 houses. Each homestead has domestic 
animals including cows, goats, sheep, chickens, and dogs. Ani- 
mal pens and corrals for animals are built adj acent to or near 
the buildings in each homestead, and in some cases the build- 
ings were within the corrai. The economie status of the popu- 
lation of Parkarin is uniformly low. Most villagers live in huts 
made of sticks and mud, tend to domestic animals, and have 
no income. Terrain throughout the village is very dry with 
scattered Acacia trees, thorn bushes, and little ground veg- 
etation. The population of Loboi is a mix of very poor fami- 
lies, day laborers, shopkeepers, and professional staff working 
at the Lake Borgoria National Park. Loboi has a health post, 
a number of shops, schools, and public utihties. Parkarin has 
no health post, therefore, residents must travel 7 km to Mari- 
gat when in need of medicai assistance. These health posts are 
staffed by nurses and community health workers who are 
responsible for patient screening for kala azar, as well as tak- 
ing blood smears and provi ding malaria treatment to patients. 
Although the workers are famihar with kala azar cUnically, 
there are no tools for diagnosing or treating this disease in the 
endemie area. Another dramatic difference between the two 
villages is the amount of vegetation. Loboi has more over- 
head cover from Acacia trees, whereas Parkarin has more 
open areas with denser scrub and a tree Une of Acacia con- 
centrated along the irrigation canal. Loboi also has more ter- 
mite mounds (A^ = 188; 87.1 mounds/km^) than Parkarin {N 
= 113; 33.7 mounds/km^). 

Study population. This study included 489 (53%) research 
volunteers from the population of Loboi and Parkarin. It was 
composed of individuals of both sexes, 5 years of age or older, 
that volunteered for the study. Ali participants had been liv- 
ing in the community for at least 3 months and called the 
community their primary home. This study was approved by 
the Walter Reed Army Institute of Research (WRAIR) Hu- 



man Use Review Committee (HURC) and the Kenya Medi- 
cai Research Institute (KEMRI) Ethical Review Committee. 
Informed written consent was obtained from ali volunteers 
prior to the employment of any examining, sampling, and 
questioning measures. Written parental consent was obtained 
for children under 18 years old. Children under the age of 5 
years and temporary visitors to the communities were ex- 
cluded from the study. 

Survey components. A community-wide survey consisting 
of a census, GIS mapping, blood sampling by venipuncture, 
medicai histories, and epidemiologie questionnaires was con- 
ducted in the two villages during the period from May to July 
2001. The epidemiologie questionnaire captured Information 
concerning house construction; the presence and location of 
domestic animals; the presence of domestic dogs; the prox- 
imity of the resident from typical sand fly breeding sites (ter- 
mite mounds and animai enclosures); how long the family had 
lived at the present site; and, the daytime occupation of each 
resident. Each participant was assigned a unique personal 
identifier and household address. Ali participant Information 
was recorded in a database and linked to diagnostics results. 

Medicai examination. A team composed of two clinicians 
experienced in VL diagnosis, a guide, and several clinical as- 
sistants and translators worked systematically through each 
village. Whereas the clinician conducted the medicai exami- 
nation, an assistant administered the medicai history ques- 
tionnaire. The study clinician recorded the participant's con- 
dition and any clinical abnormalities and/or VL-related symp- 
toms. Degree of splenomegaly was determined by palpitation 
and graded by use of the Hackett model. ^' 

Serology. Seroprevalence was determined by an ELISA 
based on an antigen derived from Leishmania donovani 
promastigotes cultivated in a protein-free, serum-free 
medium. ^^■'^^ Leishmania donovani (Ld, Strain WR0130E) 
washed promastigotes were inoculated into 200 mL of a de- 
fined, conditioned protein-free medium (XOM) to give a final 
density of 1 x 10^ cells/mL. The parasites were incubated at 
25°C for 72 hours in roller bottles. Thereafter, the spent me- 
dium was harvested by centrifugation at 11,000 x g for 20 
minutes and the relative protein concentration of the soluble 
antigens was estimated by measuring the optical density at 
280 nm.2° 

Piate sensitization was effected by coating polystyrene, 96- 
well microtiter plates (Immulon 4, Dynatech Laboratories, 
Chantilly, VA) with 100 |xL of the Ld antigen solution (5 jjig 
protein per well). 

Plates were then blocked with 0.5% casein (Sigma Chemi- 
cal Co., St. Louis, MO) in PBS for 1 hour at room tempera- 
ture. The blocking buffer was removed by aspiration and the 
serum samples (100 \lL of 1:1000 dilution) and appropriate 
Controls were added to the microtiter piate and the piate 
incubated at 26°C for 40 minutes. After washing with 0.05% 
PBS-Tween-20 (PBS-Tween) buffer four times, goat anti- 
human IgG (whole molecule) conjugated with horseradish 
peroxidase (Kirkegaard & Perry Laboratories Inc., Gaithers- 
burg, MD) was added at a 1:5,000 dilution in blocking buffer 
and the piate incubated at 26°C for 30 minutes. The piate was 
then washed four times with PBS-Tween buffer and 100 jjlL of 
TMB substrate (Kirkegaard & Perry Laboratories Inc.) was 
added to each well. The piate was incubated in the dark, and 
the optical density (OD) was periodically read at 650-nm 
wavelength in an ELISA piate reader (Molecular Devices, 
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Menlo Park, CA) until the OD value of a reference positive 
control reached 0.75. At this point, 100 ^jlL of a stop solution 
(1.0 M phosphoric acid) was applied to the piate and the final 
OD reading taken at 450 nm. A positive reference serum was 
used in ali plates, and only interassay variation of less than 
10% was accepted. The lower hmit of positivity (cutoff) was 
determined by the mean of the negative controls subset + 3 
standard deviations.^^ 

Global Information System integration. Global Informa- 
tion System (GIS) and spatial clustering statistical techniques 
were applied to evaluate the presence of high-risk areas and 
spatial clustering Leishmania seroprevalence. Ali residences, 
streets, public buildings, latrines, termite mounds and other 
features of interest in the villages were mapped using Trimble 
ProXR GPS receivers to permit the calculations of latitude, 
longitude, and altitude. During the mapping, each house was 
assigned a unique household Identification number and Infor- 
mation about construction type was also recorded. Pathfinder 
Office software, version 1.1 (Trimble Navigation, Sunnyvale, 
CA) was used to perform differential correction of ali feature 
locations and to create a locational database for use in GIS 
analyses. Seroprevalences were calculated by household. For 
the GIS analysis, ArcView software, version 3.0, (Environ- 
mental Systems Research Institute, Inc., Redlands, CA) was 
used to merge the locational database with the study group 
data. A surface analysis of cumulative Leishmania risk across 
the entire village at the household level was performed by 
using the Inverse Distance Weighted (IDW) interpolator 
modeP^ using ArcView Spatial Analyst v.1.0 (ESRI, Inc. 
Redlands, CA). IDW weights the contribution of each input 
(control) point by a normalized inverse of the distance from 
the control point to the interpolated point. IDW assumes that 
each input point (household) has a locai influence that di- 
minishes with distance. It weights the points closer to the 
processing points more than those farther away. A specified 
number of points (or ali points) within a specified radius are 
used to determine the output value for each location. 

Statistical analysis. Chi-square and Fisher's exact tests were 
used to compare differences in proportions. The Student t test 
or Mann-Whitney U, a nonparametric test, was used to com- 
pare differences in continuous variables. To evaluate ELISA 
results and splenomegaly size on a qualitative scale, the Bar- 
tholomew test was applied.^^ B ecause the study population 
included sets of familial members, the data were presumed to 
violate the standard logistic regression assumption of inde- 
pendent response probabili ties across observations. Potential 
risk factors and VL seropositivity were evaluated by odds 
ratios (OR) using random effects logistic regression for each 
village, where the homesteads were defined as the group vari- 
able.^"* To describe the mount of aggregation existing in VL 
seropositivity within homesteads and/or households units, the 
percentage of explained variance attributed (rho) was esti- 
mated after adjusting for age and gender. Ali reported con- 
fidence intervals (CI) are 95%, and ali reported P values are 
two-sided. These statistical analyses were carried out using 
Stata V. 6.0 (Stata Corporation, College Station, TX). 

Spatial analysis. To evaluate the presence of spatial disease 
clusters and to identify focalized areas of "high risk" for VL, 
we applied a spatial scan statistic.^^ The spatial scan statistic 
is well suited for geographic disease surveillance. This method 
takes into account the uneven spatial geographic distribution 
of cases and population densities; it does not require a priori 



assumptions about the number, place, or size of locations that 
may be identified as clusters; it adjusts for multiple testing 
inherent in the search for multiple clusters; and it searches for 
either high or low incidence or prevalence areas. The spatial 
scan statistic works by aggregating together the unique com- 
binations of small-area geographies that have a high probabil- 
ity of being clusters. This statistical test can detect spatial 
clusters of any size located anywhere in the study area. The 
spatial scan statistic imposes a geographic circular window to 
perform purely spatial analysis of varying size on the map 
surface and allows its center to move so that at any position 
and size across the study area, the circular window includes 
different sets of adj acent households. As the circular window 
is placed at each household, this spatial method creates a 
large number of distinct geographical circles, with different 
sets of households areas within them, and each is a possible 
candidate for containing a spatial cluster of prevalent leish- 
maniasis cases. For any given geographic position of the 
household, the radius of the circular window varied continu- 
ously from zero to a previously user-defined maximum. Al- 
though the choice of maximum cluster window size is some- 
what arbitrary, it is important to make the choice of maximum 
cluster window a priori to avoid the problems of multiple 
hypothesis testing. In this analysis, we chose 50% of the total 
study population as the maximum to consider. By choosing 
50% as the maximum, we evaluated ali sizes from zero per- 
cent up to 50% and adjusting for the multiple testing related 
to each of them. We assumed that the geographic spatial 
distribution of leishmaniasis cases foUow a Bernoulli distri- 
bution. The most likely spatial cluster was determined by 
computing maximum hkelihood ratios. The spatial scan sta- 
tistic uses the Monte Carlo simulation to evaluate the statis- 
tical significance of the most hkely spatial cluster. In this 
study, the simulated P value of the statistic was obtained 
through 9,999 Monte Carlo simulations, where the nuli hy- 
pothesis of no high spatial leishmaniasis clusters was rejected 
at an a level of 0.05. Spatial risk patterns for VL was ex- 
pressed as relative risk and was estimated by household as the 
observed VL seroprevalence/expected VL seroprevalence by 
using the spatial scan statistic. These spatial analyses were 
carried out by using the SaTScan v.2.1 software.^*' 

RESULTS 

Study area and population. The census conducted at the 
onset of this study showed that there were a total of 643 
villagers in Loboi and 286 villagers in Parkarin living in 235 
homesteads: 176 in Loboi and 59 in Parkarin. Population den- 
sity in Loboi was approximately four times that in Parkarin, 
with 279.6 and 85.4 persons per km^, respectively; 52.5% of 
the villagers were males, and 62.2% were less than 20 years 
old with a median age of 15 years. The mean number of 
persons per homestead in Loboi was 3.7 (standard deviation 
[SD] = 2.7) and in Parkarin it was 4.8 (SD = 3.8). There was 
no significant difference in mean age and gender between 
both villages. In wall construction, there was a clear and sig- 
nificant difference in both villages. Parkarin had more houses 
made of mud and sticks than Loboi (78.0% versus 31.4%, P < 
0.05) (Table 1). Residents of Parkarin had lived an average of 
11.5 years (SD = 8.9) in their community, whereas residents 
of Loboi had lived in theirs an average of 9.8 years (SD = 
10.7). People generally resided longer in their houses in 
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Table 1 
Demographic characteristics of ali population and participants tested by study sites in Kenya, 2001 







I 


i'opulation 






Fested 


Feature 


AH people No. 


Loboi No.(%) 


Parkarin No. (%) 


Ali tested No. (%) 


Lobi No. (%) 


Parkarin No. (%) 


Population/homestead 














No. of inhabitants 


929 


643 (69.2) 


286 (30.8) 


490 (52.7) 


290 (59.2) 


200 (40.8) 


No. of homesteads 


235 


176 (74.9) 


59 (25.1) 


180 (76.6) 


127 (70.6) 


53 (29.4) 


Wall construction (stick and mud) 


101 


55 (31.4) 


46 (78.0)* 


88 (87.1) 


45 (35.7) 


43(81.1)* 


Mean (SD) number of people per homestead 


4.0 (3.0) 


3.7 (2.7) 


4.8 (3.8) 


2.7 (2.1) 


2.3(1.7) 


3.8 (2.7)* 


Gender 














Female 


439 


305 (47.7) 


134 (46.9) 


245 (55.8) 


154(53.1) 


91 (45.5) 


Male 


486 


334 (52.3) 


152 (53.1) 


245 (50.4) 


136 (46.9) 


109 (54.5) 


Age group (years) 














< 10 


287 


201 (33.2) 


86 (32.8)* 


81 (28.2) 


47 (17.2) 


34 (18.4)* 


10-20 


252 


167 (27.6) 


85 (32.4) 


176 (69.8) 


99(36.1) 


77 (41.6) 


21-30 


184 


144 (23.8) 


40 (15.3) 


106 (57.6) 


74 (27.0) 


32 (17.3) 


31-40 


87 


63 (10.4) 


24 (9.2) 


58 (66.7) 


38(13.9) 


20 (10.8) 


>40 


57 


30 (5.0) 


27 (10.3) 


38 (66.7) 


16 (5.8) 


22(11.9) 


Mean (SD) age (years) 


18.0(14.1) 


17.9 (13.6) 


18.4(15.2) 


21.3(12.7) 


21.4(12.3) 


21.3(13.4) 



Note: Denominators total varied slightly due to missing data; SD, standard deviation. 
^' P < 0.05 by chi-square or Fisher's exact test compared between villages. 



Parkarin (mean 4.5 years; SD = 4.5) than in Loboi (mean 2.8 
years; SD = 3.3). 

Sample population. From the potential study population, 
289 villagers in Loboi and 200 in Parkarin volunteered to 
participate in this study, representing 44.9% and 69.9%, 
repectively. Of these, 245 (50%) were males (Table 1). The 
participants Uved in 180 homesteads (76.6% of the total 235 
homesteads). The population studied in Loboi was not sig- 
nificantly older than that of Parkarin (P = 0.924). The aver- 
age number of people per homestead sampled in Loboi was 
significantly lower than in Parkarin {P < 0.001). There were 
also significant differences with respect to wall construction, 
domestic animals, presence of goats, dogs and where animals 
were kept at night between villages (P < 0.05) (Table 2). 

Clinical and serological findings. Of the 481 (98.4%) indi- 
viduals examined, a total of 155 (32%) presented with some 
degree of splenomegaly: 90 (45% of the 200 participants) in 
Parkarin and 65 (23%) in Loboi. Of these 155 individuals, 
only 11 (7.1%) presented with other classic signs and symp- 
toms of kala azar such as pallor, weight loss, and fever. How- 
ever, we were unable to rule out other potential causes of 
chronic fever such as typhoid fever or malaria. 

The overall antibody prevalence by ELISA was 31.5% (154 
of 489). The antibody prevalences were significantly different 
between the two villages, 52.5% and 16.9% in Parkarin and 
Loboi, respectively (P < 0.001). The difference in seropreva- 
lence between females (29.4%, 72 of 245) and males (33.6%, 
82 of 244) was not statistically different (P = 0.364), even 
when the analysis was performed for each village separately. 
The mean age of positive subjects (25.6 years) was greater 
than that of subjects without antileishmanial antibodies (19.6 
years). The difference was statistically significant whether the 
population was analyzed as a whole (P < 0.001) or separately 
for each village. There was a strong linear association of in- 
creasing seroprevalence of VL with age groups (P < 0.001 by 
X^ for trend) in both villages (Figure 1). An association (P = 
0.022 by Bartholomew test) between ELISA results and spleno- 
megaly size on a qualitative scale was also found (Table 3). 

In Loboi, a significant level of aggregation was found only 
within homesteads units (rho = 27%, 95% CI = 13 to 48, P 
< 0.001) and not within household units. Households or 



household aggregation did not explain a significant propor- 
tion of the variance of VL seropositivity in Parkarin. 

Epidemiologìe faetors. From the serological and epidemio- 
logie surveys, we found a number of faetors strongly associ- 
ated with altered risk of leishmaniasis. For each 10-year in- 
crease in age, the likelihood of being seropositive was higher 



Table 2 
Village comparison using certain epidemiological faetors of home- 
stead by study sites in Kenya, 2001 



Feature 


Visited No. 


Loboi No. (%) 


Parkarin No. (%) 


No. of homesteads 








in the village 


180 


127 (70.6) 


53 (29.4) 


Presence of domestic animals 






No 


71 


56 (45.5) 


15 (28.3)* 


Yes 


105 


67 (54.5) 


38 (71.7) 


Presence of cattle 








No 


114 


81 (65.9) 


33 (62.3) 


Yes 


62 


42 (34.1) 


20 (37.7) 


Mean (SD) number 








of cattle 


2.5 (4.9) 


2.8 (5.6) 


1.9(3.1) 


Presence of sheep 








No 


156 


112(91.1) 


44 (83.0) 


Yes 


20 


11 (8.9) 


9 (17.0) 


Mean (SD) number 








of sheep 


1.4 (8.3) 


1.7 (9.8) 


1.2 (2.9) 


Presence of goats 








No 


99 


79 (64.2) 


20 (37.7)* 


Yes 


77 


44 (35.8) 


33 (62.3) 


Mean (SD) number 








of goats 


7.6 (12.4) 


5.7 (17.3) 


12.2 (12.4)* 


Presence of chicken 








No 


91 


66 (53.7) 


25 (47.2) 


Yes 


85 


57 (46.3) 


28 (52.8) 


Mean (SD) number 








of chicken 


3.5 (5.2) 


3.6 (5.5) 


3.3 (4.5) 


Presence of dogs 








No 


127 


96 (78.7) 


31 (58.5)* 


Yes 


48 


26 (21.3) 


22 (41.5) 


Animai shelters 








Inside house 


23 


21 (17.4) 


2 (4.0)* 


Adjacent corrai 


76 


44 (36.4) 


32 (64.0) 


Away from homestead 


2 


1 (0.8) 


1 (2.0) 


Note: Denominators total varied slightly due to 


missing data: SD, t 


ìtandard deviation. 


=" P < 0.05 by chi-square or Fisher 


's exact test compared between villages. 
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Figure 1. Seroprevalence of Leishmania-speciììc antibody by age 
group and study sites in Kenya, 2001. 



in Loboi (OR = 1.88, 95% CI = 1.32 to 2.69, P < 0.0001) than 
Parkarin (OR = 1.52, 95% CI = 1.16 to 1.99, P = 0.002). In 
univariate logistic regression analysis for Loboi, residents liv- 
ing in households made of sticks and mud, as well as presence 
of domestic animais, especially cattle, sheep, and goats, were 
more likely to be seropositive to VL. In the multiple logistic 
regression analysis adjusted for age and gender, these poten- 
tial risk factors remained statistically significant associated to 
VL with the exception of cattle. In Parkarin, presence of do- 
mestic animais in households (especially sheep and goats) was 
negatively associated with a positive ELISA and appears to 
be protective in the univariate analysis (Table 5). However, 
only presence of goats remained negatively associated to VL 
in the multiple logistic regression analysis. Therefore, these 
data suggest an inverse odd ratio risk relationship in the two 
villages studied. In both villages gender, proximity to active or 
inactive termite mounds, and presence of dogs were not sig- 
nificantly associated with VL. 

Spatìal analysis. Based on the serological findings, we gen- 
erated risk pattern maps by household for Loboi (Figure 2) 
and Parkarin (Figure 3) that defined areas of high risk. The 



Table 3 
Splenomegaly size of visceral leishmaniasis by ELISA 







ELISA 






Spleen size category 


Positive no. 




Negative no. 


Total no. 


None 


89* 




237 


326 


Slight 


30 




38 


68 


Palpable 


28 




48 


76 


Considerable 


3 




5 


8 


Severe 


2 




1 


3 


Total 


152 




329 


481 



ELISA, enzyme-linlied immunosorbent assay. 

* /• = 0.022 by Bartholomew test (association between ELISA results and splenomegaly 
size on a qualitative scale}. 



spatial scan statistic revealed the presence of a spallai cluster 
of Leishmania-specific seropositives in the village of Loboi 
(Figure 4). This spatial cluster contained 11% (16 of 152) of 
the total houses sampled in the village and 29% (14 of 49) of 
the total ELISA cases. Villagers living within this spatial clus- 
ter were 3.3 times more likely to have Leishmania-specific 
antibodies than people living in other areas in the village 
(Table 4; P = 0.0003). In addition, this spatial cluster 
occurred in the northwestern corner of the village in a 
wooded area in very dose proximity to a seasonal river, a 
kind of setting commonly associated with VL. The resi of 
the village consisted of more open areas with scattered 
Acacia trees, shrubs, and little ground cover. There were two 
other spatial clusters that appear graphically (bui not statis- 
tically) in Loboi, and their settings were nearly identical 
to that of the primary cluster. There was not a specific high 
risk pattern and no significant spatial cluster was found 
(P = 0.296). The spatial distribution of seropositives among 
households appeared to be homogeneous (Figure 3), and this 
may be due to the high seroprevalence of VL found in this 
village. 

DISCUSSION 

To our knowledge, this study represents the first geo- 
graphic epidemiologie analysis in two villages to study the 
microepidemiology of VL in Africa and probably in other 
continents. We identified microgeographic areas where there 
is more likely to be high Leishmania exposure and spatial 
cluster of seropositives through GIS mapping and spatial sta- 
tistic techniques. 

Baringo Districi has been the focus of many epidemiologie 
studies on leishmaniasis. As early as 1983, Jahn and Dies- 
feld^^ used a crude ELISA for detecting Leishmania specific 
antibodies in VL patients seen in the Baringo Districi Hos- 
pital al Kabarnet. With limited success, they applied this 
ELISA in clinical routine and sero-epidemiologic surveys. 
Low titers were recorded in healthy individuals from VL foci, 
but their values were easily distinguished from those of active 
patients. After that small study, Jahn and co-workers, using 
the same ELISA, showed that children between 2 and 15 
years old were the most affected age group for VL and that 
male patients predominated shghtly al 57%.^^ Ali VL cases 
came from the semiarid and arid parts of the districi below 
1500 m, where pastoralism predominated. A positive corre- 
lalion was reported for active cases and the proximity of 
homes to seasonal rivers. However, no satisfying explanation 
was found for the clustering of cases. ELISA values above the 
cutoff, taken as the borderline nonspecific reaction, were 
found in about half of the population, suggesting that asymp- 
tomatic infections were common.^'' 

Schaefer and coUeagues tested Baringo District inhabitants 
from 26 clusters, averaging 97 people, with the leishmanin 
skin test (LST).^° There was an obvious aggregation of LST 
positivity centered on recent VL cases with approximately 
11% of the 2,411 individuals tested having a marked DTH 
reaction. The level of LST positivity was twice as high in 
males as females. LST positivity increased with age to a stable 
level after 15 years of age, reflective of an endemie situation. 
Schaefer and co-workers later measured Leishmania-specific 
antibodies using the direct agglutination test (DAT) in a de- 
fined, endemie rural area of Baringo District composed of 30 
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Figure 2. Spatial risk patterns for visceral leishmaniasis in Loboi, Kenya, 2001. This figure appears in color at www.ajtmh.org. 



clusters, each averaging 98 people.^^ From 2,934 individuals 
samples, 78 were DAT seropositive; 54 of those had a history 
of VL. The seroconversion rate was 9 of 1,000 person-years of 
observation among 2,332 seronegative individuals retested 
the following year. During the entire study period, VL was 
diagnosed in only 10 patients, with an incidence rate of 2.2 of 
1,000 person-years of observation. Household contacts of in- 
dividuals with previously confirmed VL had a higher fre- 
quency of DAT positivity than the rest of the population. 
These results suggested domiciliary transmission. 

Using another approach, Schaefer and colleagues used 
PCR with capillary blood samples dried on filter paper.^^ 
Assaying 20 parasitologically confirmed VL case-patients, 21 
subclinical cases, and 11 healthy controls from a longitudinal 
study of anthroponotic VL in B aringo District, they were able 
to detect Leishmania DNA 10 months before diagnosis and 
up to 3 years after treatment that was classified as successful. 
Alarmingly, these findings showed that subclinical and 
treated cases may remain potential reservoirs for long peri- 
ods. 

The study results suggest that males and females have simi- 
lar risks to VL, and seropositive cases are clustered within 
households. The results of this study are interesting and, at 
the same time, confounding. We see dramatic differences in 
seroprevalences between the two villages (16.9% versus 



52.5%). GIS mapping and spatial scan statistic showed the 
presence of a spatial cluster only in the village with the lowest 
seroprevalence. This spatial cluster found may be related to 
the wooded environment of the northwestern corner of Loboi 
that provides a good habitat for sand flies.^'' When merging 
ELISA data with the results of the epidemiologie survey, we 
find strong positive associations between seropositivity, age, 
and poor living conditions. House building materials (mud 
and sticks) related to lower socioeconomic status were asso- 
ciated with higher risk. However, we find that the presence of 
domestic animals is only associated with increased risk in 
Loboi. The figures from Parkarin indicate that there was an 
inverse relationship between the presence of animals and ex- 
posure to Leishmania. For both villages, there was a signifi- 
cant association between ELISA seroprevalence and distance 
to the nearest corrai (P = 0.015 for Loboi and P = 0.004 for 
Parkarin) (data not shown). 

We surmise that this result may be due to the fundamental 
differences between the two villages with respect to proximity 
of rivers, vegetation cover, population density, and economie 
make-up of the population. However, there may be limita- 
tions in the statistical analysis due to the extremely high se- 
roprevalence found in Parkarin and the small size of this vil- 
lage. Our data suggests strategies for targeted control and for 
prioritization of scarce resources. A community-based VL 
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Figure 3. Spatial risk patterns for visceral leishmaniasis in Parkarin, Kenya, 2001. This figure appears in color at www.ajtmti.org. 
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Figure 4. Map showing location of the most likely spatial cluster of visceral leishmanisis in Loboi. Kenya, 2001. This figure appears in color 
at www.ajtmh.org. 
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Table 4 
Most likely spatial cluster of visceral leishmaniasis by study sites in 
Kenya, 2001 



Feature 


Loboi 


Parkarin 


Study area 
No. of houses 


152 


75 


No. of inhabitants 


290 


200 


No. of positive inhabitants by ELISA 
Seroprevalence 
Most likely spatial cluster 
No. of houses 


49 
16.9 

16 


105 
52.5 

5 


No. of inhabitants 


25 


11 


No. of positive inhabitants by ELISA 

Seroprevalence 

Overall relative risk* 


14 
56.0 
3.28 


10 

90.9 
1.73 


P value 


0.0003 


0.2961 



ELISA, enzyme-linked immunosorbent assay. Seroprevalence was estimateci as the num- 
ber of positive inhabitants by ELlSA/total inhabitants tested. 

^' Relative risk = observed VL seroprevalence/expected VL seroprevalence. 



control or suppression program could be formulated on edu- 
cating residents of this endemie area about the risk associated 
with house construction and the proximity of domestic ani- 
mais to one's living quarters. These points are likely related to 
sand fly habitat, whereby the sand flies seek out a favorable 
habitat in this arid setting that would be supportive of egg 
laying and survival of immature forms.^'' Most residents kept 
their animals at night in corrals no more than 30 m from their 
homes. A higher and statistically significant difference (P = 
0.001) of the ratio of animals to people was found in Parkarin 
(median = 17.5) than in Loboi (median = 5). The animals 
serve as a blood source,^^"^'^ and the accumulation of animai 
dung may be attractive to the sand fly, drawing the vectors 
into closer association with the humans and increasing the 
risk of being bitten. In addition, the edges of the corrai con- 
structed from cut thorn bushes provide good areas for rodent 
burrows and nesting sites. Positive correlation of disease and 
the presence of domestic animals have been shown in some 
studies,^^ whereas others have shown no significant effect^^ or 
have noted an inverse relationship between domestic animals 
and disease presence. ^^ However, the nature of the relation- 
ship between disease transmission and domestic animals is 
complex."^'^^ Among the factors that can affect feeding be- 
havior, host odor, heat loss, and CO2 production are impor- 
tant stimuli for orientation of blood sucking insects in the 
open (i.e., nonforest) situations,^^ diverting blood feeding 
from people to animals. Blood meal studies of P. martini, the 
vector of kala azar in the study area, has shown it to have a 
broad host range among domestic animals, with a preference 



for goats.^"* Zooprophylaxis explains the situation in Parkarin 
where the higher ratio of goats appears to have a protective 
effect on the human population. 

VL is a grave public health problem in this area that im- 
poses an additional strain on the locai health authorities and 
is unlikely to be resolved by the current strategies. Under- 
standing the transmission dynamics of VL could lead to sus- 
tainable prevention and control measures. Alternative meth- 
ods of vector control, other than the conventional indoor 
spraying of houses with residuai insecticide (which can be 
prohibitively expensive) should be considered. Repeated pes- 
ticide apphcations to the whole community can reduce sand 
fly numbers, but as the reduction is short hved, this method is 
used only in epidemics."*^ 

The termite hiils associated with leishmaniasis transmission 
in Kenya^® are common throughout the two villages. Ai- 
though no association was shown in this study, the termite hili 
habitat is the favored breeding and resting sites of P. mar- 
tini.'^'^''^^ This points to targeted control strategies such as in- 
secticidal apphcations to resting habitats such as termite 
mounds'*'^ and insecticide barrier spraying.'*^ The feeding pro- 
pensity of P. martini on livestock and the efficacy of zoopro- 
phylaxis could be enhanced by using the livestock not only to 
divert the vector from humans to animals but to attract them 
to contact with insecticide-treated livestock.'*^''*'' Because P. 
martini is active during the night when people sieep, the use 
of treated bed nets to block transmission couid provide con- 
siderable protection.''' A prospective intervention study 
would be needed to evaluate the effectiveness of targeting con- 
trol interventions at high-risk areas identified by tiiis study. 
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Table 5 
Potential risk factors associated with visceral leishmaniasis seropositivity by study sites in Kenya, 2001 







L 


3boi 






Parkarin 




Feature 


OR (95% CI) 


P value 


AOR (95% CI) 


P value 


OR (95% CI) 


P value 


AOR (95% CI) 


P value 


Wall construction type 


















Stone, wood. sheet metal (sticks 


















only, mud) 


4.3 (2.1-8.9) 


< 0.001 


5.5 (2.4-12.5) 


< 0.001 


2.0 (0.9-4.2) 


0.083 


1.6 (0.8-3.2) 


0.165 


Presence of domestic animals (no.) 


3.8 (1.8-8.2) 


0.001 


4.1 (1.6-10.4) 


0.003 


0.2 (0.1-0.5) 


0.001 


0.4 (0.2-0.8) 


0.016 


Presence of cattle (no.) 


2.0 (1.1-3.7) 


0.036 


1.9 (0.9-4.6) 


0.103 


0.6 (0.3-1.1) 


0.091 


0.7 (0.4-1.4) 


0.387 


Presence of sheep (no.) 


2.6 (1.2-5.7) 


0.014 


3.3 (1.2-8.9) 


0.018 


0.4 (0.2-0.9) 


0.026 


0.6 (0.3-1.3) 


0.178 


Presence of goats (no.) 


3.5 (1.8-7.0) 


< 0.001 


4.5 (2.1-9.9) 


< 0.001 


0.3 (0.2-0.7) 


0.003 


0.4 (0.2-0.8) 


0.016 



ELISA, enzyme-linked immunosorbent assay; OR, odds ratio; AOR, adjusted odds ratio by gender and age; 95% CI, 95 percent confidence interval. Categories in parentheses describe the 
reference group for odds calculations. 
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